Minimization of phonon-tunneling dissipation in mechanical resonators 
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Micro- and nanoscale mechanical resonators have recently emerged as ubiquitous devices for use 
in advanced technological applications, for example in mobile communications and inertial sensors, 
and as novel tools for fundamental scientific endeavors. Their performance is in many cases limited 
by the deleterious effects of mechanical damping. Here, we report a significant advancement towards 
understanding and controlling support-induced losses in generic mechanical resonators. We begin 
by introducing an efficient numerical solver, based on the "phonon-tunneling" approach, capable 
of predicting the design-limited damping of high-quality mechanical resonators. Further, through 
careful device engineering, we isolate support-induced losses and perform the first rigorous experi- 
mental test of the strong geometric dependence of this loss mechanism. Our results are in excellent 
agreement with theory, demonstrating the predictive power of our approach. In combination with 
recent progress on complementary dissipation mechanisms, our phonon-tunneling solver represents 
a major step towards accurate prediction of the mechanical quality factor. 



Mechanical coupling of a suspended structure to its 
supports is a fundamental energy loss mechanism in 
micro- and nanomechanical resonatorgl". Referred to 
variously as clamping^ or anchor losa^, this process re- 
mains significant even in devices fabricated from high- 
quality materials operated in vacuum and at cryogenic 
temperatures, and is in fact unavoidable in any non- 
levitating system. Although much progress has been 
made towards the understanding of mechanical dissipa- 
tion at the micro- and nanoscald^li]^ obtaining reliable 
predictions for the fundamental design-limited quality 
factor Q remains a major challenge while direct exper- 
imental tests are scarce. At the same time, the im- 
plementation of high-quality micro- and nanomechanical 
systems is becoming increasingly important for numer- 
ous advanced technological applications in sensing and 
me trol ogy, with select examples including wireless fil- 
tera^, on-chip cl ocks^ , microscop}EHini and molecular- 
scale mass sensing—'--', and recently for a new genera- 
tion of macroscopic quantum experiments that involve 
mesoscopic mechanical structure4i2tt20|_ jjere we intro- 
duce a finite-element enabled numerical solver for cal- 
culating the support-induced losses of a broad range of 
low-loss mechanical resonators. We demonstrate the effi- 
cacy of this approach via comparison with experimental 
results from microfabricated devices engineered in order 
to isolate support-induced losses by allowing for a sig- 
nificant variation in geometry, while keeping other res- 
onator characteristics approximately constant. The effi- 
ciency of our solver results from the use of a perturbative 
scheme that exploits the smallness of the contact area, 
specifically the recently introduced "phonon-tunneling" 
approactpH. This results in a significant simplification 
over previous approaches and paves the way for CAD- 
based predictive design of low-loss mechanical resonators. 

The origins of mechanical damping in micro- and 
nanoscale systems have been the subject of numerous 
studies during the last decades and several relevant 
mechanisms for the decay of acoustic mechanical exci- 



tations, i.e. phonons, have been investigatecpl^. These 
include: (i) fundamental anha rmonic effects such as 
phono n-ph onon interaction d^^^^l , thermoelastic damping 
(TEDJ^IMISII^ and the Akhiezer effedP^l (ji) viscous 
or fluidic damping, involving interactions with the sur- 
rounding atmosphere or the compression of thin fluidic 
layera^SHlS Qfi) materials losses driven by the relaxation 
of intrinsic or extrinsic defects in the bulk or surface of 
the rcsonatop22n23l for which the most commonly stud- 
ied model is an environment of two-level fluctuatorP^M] 
and (iv) support-induced losses, i.e. the dissipation in- 
duced by th e una voidable coupling of the resonator to 
the substrate^'^^SMSil^ which corr espon ds to the radiation 
of clastic waves into the supporta^323HIl, This last mech- 
anism poses a fundamental limit as vibrations of the sub- 
strate will always be present. 

These various dissipation processes add incoherently 
such that the reciprocals of the corresponding Q-values 
satisfy 1/Qtot = J2i ^/Qi^ where i labels the different 
mechanisms. Thus, in a realistic setting care must be 
taken to isolate the contribution under scrutiny. In con- 
trast to all other damping mechanisms (i)-(iii) which ex- 
hibit various dependencies with external physical vari- 
ables such as pressure and temperature, support-induced 
dissipation is a temperature and scale-independent phe- 
nomenon with a strong geometric character that is 
present in any suspended structure. Moreover, its scale- 
independence implies that the same analysis can be ap- 
plied to both micro- and nanoscale devices. We exploit 
this geometric character in order to isolate the support- 
induced contribution and obtain a direct experimental 
test of phonon-tunneling dissipation. 

We first introduce our numerical solver, which provides 
a new technique to efficiently model support-induced 
losses for a broad class of mechanical structures. Pre- 
vious approaches have relied on either the direct solution 
of an el astic wave radiation problem involving the sub- 
strate^SnSIlSl or t he sim ulation of a perfectly absorbing 
artificial boundarySHlll^ with applications typically lim- 



ited to a few specific geometries^ZEIIllllSl. Jn contrast, 
our technique represents a substantial simplification in 
that it reduces the problem to the calculation of a per- 
fectly decoupled resonator mode together with free elas- 
tic wave propagation through the substrate in the ab- 
sence of the suspended structure. A key feature of our 
method is to combine a standard finite-element method 
(FEM) calculation of the resonator mode together with 
the use of an extended contact at the support. This al- 
lows us to treat complex geometries taking proper ac- 
count of interference effects between the radiated waves. 

In analogy to radiation tunneling in photonics and elec- 
tron tunneling in low-dimensional structures, we adopt 
a "phonon tunneling" picture to describe the support- 
induced losses^B. In this picture the mechanical reso- 
nance of interest, characterized by frequency ujr^ is re- 
garded as a phonon cavity that is weakly coupled to the 
exterior by a hopping process, whereby the elastic energy 
leaks out of the resonator through the narrow contact 
areas from which it is suspended. Within this frame- 
work, one can start from the harmonic Hamiltonian asso- 
ciated with the elastic scattering eigenmodes of the entire 
structure, including the substrate, and derive a quantum 
model for the Brownian motion experienced by each res- 
onance of the suspended structure. 

The corresponding weak tunnel couplings can be ob- 
tained to lowest order in the small parameter k^d, where 
1/fcfl is the characteristic length scale over which the res- 
onator mode varies appreciably and d is the characteristic 
dimension of the contact area 5* from which the resonator 
is suspended. For typical structures that exhibit high-Q 
mechanical resonances, kjid << 1 is comfortably satis- 
fied. This justifies the weak coupling approximation and 
leads to a general expression for the associated dissipa- 
tion 1/Q in terms of the overlaps between the scattering 
modes and the resonator mode. In the limit d — > the 
leading contribution is obtained by replacing the scatter- 
ing modes by the free (unperturbed) modes of the sup- 
ports, which yields^^ 
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Here <t'^ and w'^ are the stress and displacement fields 

associated with the normalized resonator mode, crq and 

Uq are the analogous fields for the continuum of support 
modes labeled by q [eigenfrequencies oj{q)], and ps and 
yOfl are, respectively, the densities of the substrate and 
resonator materials. The resonator mode should sat- 
isfy either (i) free or (ii) clamped boundary conditions 
at the contact area S depending on the behavior of the 
eigenmode when S is small, while the unperturbed sup- 
port modes should satisfy the converse. These homoge- 
neous boundary conditions correspond, respectively, to 
dS ■ (t'j^ = and m'^ = so that only one of the two 
terms in the surface integral is finite. Physical examples 
of case (i) are pedestal geometries, such as microspheres. 
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FIG. 1: Mapping out phonon-tunneling dissipation in a free- 
free resonator, (a) Schematic diagram of the resonator geom- 
etry, (b) Normalized squared center of mass displacement of 
a single auxiliary-beam-central-resonator contact calculated 
via FEM (the inset shows the profile of the free-free mode 
as approximated by Euler-Bernoulli theory), (c) Simulated 
dissipation [cf. Eq. [2] as a function of the auxiliary beam's 
j/-coordinate (ya)- Values corresponding to 8 discrete ge- 
ometries were calculated here with t = 6.67 /im, Ws = 7 pra, 
w = 42 /im, L — 132 /xm, R = 116 /^m, and Lund = 27 /im — 
the line is simply a guide for the eye. The FEM calculated 
mode shapes correspond to the three extreme examples of 
the resonator design, from left to right: auxiliary beams near 
the resonator center (j/a = 13 /im), beams near the ideal 
nodal position {ya = 37.4 /xm), and beams attached at the 
ends {ya ~ 62.5 /im). The theoretical clamping loss limit 
1/Qth for nodal positioning is always finite with the geome- 
try closest to this position (indicated by the arrow) yielding 



1/Qt 



2 X 10" 



microdisks, or microtoroids, when the pedestal has good 
impedance match with the substrate and is thus included 
in the support. On the other hand examples of case (ii) 
include the planar structures investigated here, where the 
supports consist of the portion of the structure that is 
not free-standing. Further details about the scope and 
validity of Eq. (II]) are given in the SI Appendix lAJ 

Though the aforementioned framework is completely 
general, in order to investigate the predictive power of 
our approach, we focus specifically on the fiexural modes 
of a symmetric plate geometry of thickness t that is in- 
scribed in a circle of radius R, with the contact area 
S corresponding to the outer rim of an idealized cir- 
cular undercut (undercut distance of Lund)- To calcu- 
late the theoretical Q-values of such devices via Eq. ^ 
we have developed a numerical solution technique that 
determines the resonator eigenmode and eigenfrequency 



via FEM (with u^ = at 5') and uses a decomposition 
into cylindrical modes for the support, which is approx- 
imated by the substrate modelled as an isotropic elas- 
tic half-space. The latter approximation is expected to 
be quantitatively precise for the low-lying flexural reso- 
nances when the underetched gap between the suspended 
structure and the substrate satisfies h < R (where h is the 
gap height), and the largest resonant wavelength for elas- 
tic wave propagation in the substrate is smaller than the 
relevant length scales characterizing the mounting of the 
sample (see below). The aforementioned weak coupling 
condition, kfjd << 1, follows in this case from t -^ R. 
From Eq. (IT]) we obtain [cf. SI Appendix [A| for details of 
this derivation! 
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where we use cylindrical coordinates and intro- 
duce the dimensionless displacements Un,-y{q,Vs) — 
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27rjp d0sin6'|M^,^.j,(O,i^s)pJ^(gsin0) and the hnear 
stress Fourier components /z_„ = J^dS ■ cr'j^ ■ ze™*^ with 
n = 0, ±1,±2,.... Here j = I, t, s labels the differ- 
ent types of relevant plane- wave modes u)jJy..,{f,Vs) of 
the half-space^ [i.e. longitudinal (/), transverse SV (i), 
and surface acoustic waves (s) given that transverse SH 
waves do not contribute] with Cj the corresponding speed 
of sound — as determined by the density ps , Poisson ra- 
tio i^s, and Young's modulus Eg of the substrate. We 
adopt spherical coordinates for the incident wave vector 
q with polar angle 6. It is straightforward to generalize 
the above to in-plane modes and to plate geometries in- 
scribed in a rectangle. Furthermore, we stress that the 
rim need not be continuous, as in cases where the res- 
onator volume contacts the support at a disjoint set of 
small areas. 

For experimental verification of our solver, we have 
developed "free-free" micromechanical resonators con- 
sisting of a central plate (resonator) of length L and 
width w suspended by four auxiliary beams as depicted 
in Fig. 1(a). These structures are etched from a high- 
reflectivity monocrystalline distributed Bragg reflector 
(DBR) — as described in the Methods Section, suited for 
Fabry-Perot-based optomechanical systems^. The de- 
vices used in this study constitute a variant of the previ- 
ously demonstrated free- free flexural design in which aux- 
iliary beams with widths w^ ^ w and lengths Lg = At/4 
(where A* is the resonant wavelength for the propagation 
of torsional waves) placed at the nodes of the central 
resonator mode provide noise filters to suppress support- 
induced losses^. A major drawback with the A(/4-beam 
design is that the resulting auxiliary beam length can be 
excessive. In fact for the eigenfrequencies investigated in 
this work, the corresponding beam length (> 400/im at 
1.7MHz) leads to a proliferation of low frequency flexural 
resonances that compromise the stability of the optical 
cavity and render mode identification difficult. We cir- 
cumvent this issue by utilizing instead a reduced auxil- 



iary beam length Lg ^ Af/4 chosen to avoid matching 
flexural resonances. 

The free-free design provides an ideal platform to iso- 
late and measure phonon tunneling dissipation: first, by 
altering the attachment position of the auxiliary beams, 
this design allows for a significant variation of geometry, 
while approximately preserving the frequencies and effec- 
tive surface-to-volume ratios of the resonators. As these 
characteristics are kept constant, one can rule out the 
influence of additional damping mechanisms (specifically 
those driven by internal losses and surface effects) on the 
variation in Q and hence isolate support-induced losses 
in the measured devices. Second, the free-free resonators 
provide an intuitive illustration of the strong geometric 
character of support-induced dissipation. Intuitively, the 
clamping loss will be proportional to the elastic energy 
radiated through the auxiliary beams which should ap- 
proximately scale as the squared deflection of their con- 
tacts with the central resonator [cf. Fig. 1 panels (b) and 
(c)]. Thus varying the contact position of the auxiliary 
beams results in a characteristic modulation of the damp- 
ing rate which approximately maps out the central res- 
onator mode shape. As expected, the minimum-loss de- 
sign corresponds to the geometry in which the auxiliary 
beams are attached at the nodes of the fundamental res- 
onance of the central resonator. It is interesting to note 
that the theoretical clamping loss limit 1/Qth for nodal 
positioning is always finite as described in [Fig. 1(c)]. 

For mode identification we compare the optically mea- 
sured resonator frequencies (as a function of the auxil- 
iary beam position) with the theoretical eigenfrequency 
response. The simulated values are generated using the 
geometric parameters determined via careful analysis of 
the completed resonators (cf. SI Appendix [B]). As can be 
seen in Fig. 2(b), in addition to the symmetric free-free 
resonance, there is also an antisymmetric eigenmode with 
comparable frequency. We observe no mode-coupling be- 
tween these resonances, which is consistent with the spe- 
cific mirror symmetries of the structure. The frequencies 
are accurately reproduced by the FEM simulation if we 
allow for frequency offsets that are solely dependent on 
the mode parity (262 kHz offset for the free-free mode and 
89 kHz offset for the anti-symmetric mode) . We attribute 
these shifts to a materials-related dissipation mechanism, 
involving both surface and bulk contributions (see SI Ap- 
pendix JC] for further details). 

All dissipation measurements have been performed at 
high vacuum (« lO^'' mbar) and at cryogenic temper- 
atures (20 K) in order to suppress fluidic and thermoe- 
lastic damping in the devices. Under these conditions, 
we record quality factors spanning 1.4 x 10'* to 5.1 x 10'* 
(cf. Fig. 3), with the minimum Q corresponding to the 
free-free mode of devices with an auxiliary position of 
62.5/xm and R — 116/im, and with the maximum Q to 
the geometry closest to nodal positioning (37.4^m) for 
the same radius and type of mode. For the symmet- 
ric mode, we readily observe the expected characteris- 
tic modulation in Q as a function of the placement of 
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FIG. 2: Characterization of the completed free-free resonators, (a) Optical micrograph of the 5mmx5mm chip containing the 
batch-fabricated microresonators as well as an electron micrograph highlighting a single suspended structure, (b) Simulated 
(left) and measured (right) eigenfrequencies as a function of the auxiliary beam j/-coordinate. The measured values (discrete 
points) show excellent agreement with the simulated dataset, albeit with a slight offset dependent on the parity of the mode. 
The fitting lines in the right plot correspond to a mean frequency offset of 262 kHz for the symmetric free- free modes and 
89 kHz for the neighboring antisymmetric modes (inset shows the FEM derived mode shape of the antisymmetric resonance). 
Lower panel — examples of the fitting techniques utilized for Q- value extraction including: (c) Lorentzian fitting of the free-free 
resonance (captured on a spectrum analyzer) for a device with R — 116 /im and -i/a = 29 /im resulting in Q = 4.5 x 10*; and 
(d) ringdown fitting of the same device using linear regression of the natural log of the mean-square of the free-ringdown signal 
captured single-shot with a high speed oscilloscope yielding Q = 4.46 x 10*. The inset includes the residuals to the linear fit 
showing an excellent agreement with the expected exponential decay. 



the auxiliary beams (cf. Fig. 1) with a relative varia- 
tion of AQexp/Qexp ~ 260% (~ 80%) for R = U6fim 
(R = ISl/Ltm). At the same time, the frequency variation 
is kept small "qua design" , with a range of A/// ~ 20% 
(~ 10%). In contrast the Q-values for the antisym- 
metric mode are nearly constant with Q « 2.1 x 10"' 
[Fig. 3(c)]. This is expected as the theoretical support- 
induced loss for this mode is negligible. Additionally, as 
this resonance involves mainly deformations of the aux- 
iliary beams, its dissipation is not expected to be cor- 
related with the mode shapes of the central resonator. 
The damping of this mode is instead dominated by 
other sources of dissipation, most likely by the materials- 
related losses that are also responsible for the frequency 
shifts. Thus we obtain an independent corroboration that 
the characteristic Q-variation observed for the free-free 
mode is indeed induced by the modification of the geom- 



etry rather than by the small frequency variation present 
in the devices. 

In order to quantitatively compare the measurements 
with our numerical predictions, two additional issues 
must be considered: (i) our model only captures support- 
induced losses, while other loss mechanisms may still con- 
tribute to the overall damping in the devices, and (ii) the 
parameters for the half-space model of the substrate must 
be properly chosen. Consideration (i), together with the 
fact that we have designed sets of resonators for which 
the frequencies and effective surface-to- volume ratios are 
kept approximately constant, implies that any additional 
damping mechanism that is relevant at low temperatures 
and high vacuum, but is insensitive to the variation in 
geometry, should contribute a constant offset l/Q* in 
the measured dissipation l/Qtot- Consideration (ii) is 
non-trivial given the long-wavelength nature of the elas- 
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FIG. 3: Compiled dissipation results displaying excellent 
agreement between theory and experiment, (a) and (b) Com- 
parison of experimental measurements at T = 20 K with 
theoretical dissipation values for the free-free mode of res- 
onators with measured central dimensions of 132/imx42/im 
and radius R = 116 /im and R — 131 fim, respectively. Both 
ringdown and spectrally-derived data are included, with val- 
ues averaged over two nominally-identical chips. We include 
both raw simulated data as well as fitted data (continuous 
lines are a guide to the eye) incorporating a constant offset 
1/Q* = 2.41 X 10""". For the effective substrate we utilize the 
mechanical properties of Ti, which is the main constituent 
of the positioning system on which the chips are mounted 
(p = 4540kg/m^ Es = 116 GPa, and Us = 0.34). (c) Mea- 
sured dissipation for the antisymmetric mode of the same 
structures exhibiting a lack of geometric dependence. 



tic waves radiated into the substrate. For an average 
resonator frequency « 2.12 MHz estimates of the maxi- 
mum wavelength for the freely-propagating elastic waves 
yield a value of « 2.5 mm, which largely exceeds the wafer 
thickness (300 ± 25 /im). Thus, the mechanical materi- 
als parameters for the substrate should be determined 
by the properties of the underlying stage and positioning 
mechanism in the cryostat rather than by those of the 
chip itself. Hence, we assume for the half-space the me- 
chanical properties of polycrystalline commercially-pure 
(grade 2) titanium (cf. Fig. 3), of which the bulk of the 
structure beneath the resonator consists. Taking all of 
this into account, the theory shows remarkable agree- 
ment with the measured dissipation (as shown in Fig. 3). 
It is important to note that the only free parameter used 
in the model of the free-free mode is a constant offset 
of l/Q* = 2.41 X 10^^. Although the exact nature of 
the corresponding dissipation mechanism is currently un- 
known, we assume that it arises from materials losses in 
the resonator epi-structure. As our design isolates the 
support induced-losses, these results establish the first 
systematic test of phonon tunneling dissipation in me- 
chanical resonators. 

It should be noted that most commercially viable 
resonators operate in a regime where TED dominates, 
and in some instances intuitive understandings of the 
support-induced damping^EESl have allowed for its sup- 
pression below other limiting damping mechanisms. 
Nonetheless, if current efforts to minimize TED in such 
structures at room temperature are successfulSS support- 
induced losses may pose the next challenge for max- 
imizing Q. On the other hand, in fundamental re- 
search thrusts employing high vacuum and cryogenic sys- 
tems, support-induced losses can become a dominant 
factopSZES. For example, the free-free designs explored 
here provide a route to minimize support-induced losses 
for application in optomechanical experiments utilizing 
the micromechanical resonator as an end mirror in a high- 
finesse Fabry-Perot cavity^'. To gauge the relevance of 
our "free-free" micromirror design in this context, it is 
instructive to compare the fundamental limit at nodal 
positioning Qth ~ 5 x 10^ and the maximum Q-value 
measured of 5.1 x 10**, with the corresponding results for 
the fundamental flexural mode of a clamped bridge of 
comparable dimensions. In fact for the typical dimen- 
sions considered, as required for integration in a high- 
performance Fabry-Perot cavity, we obtain a theoretical 
limit Qc~c ^ 10'^ m line with prior measurements on 
monocrystalline DBR optomechanical structureP^. 

Given the scale independent nature of support-induced 
losses, our solver applies equally well to nanoscale me- 
chanical devices. We find that for a recent demonstration 
of a nanomechanical doubly-clamped beam coupled to a 
superconducting qubit at mK temperaturea^, the mea- 
sured values for the resonator's maximum Q (« 6 x 10^) 
can be understood solely via the phonon-tunneling loss 
model (beam geometry of 0.3/xmx0.18/.tmx6/im) which 
predicts a Q of 5.4 x 10^, in excellent agreement with the 



experimental value. 

In conclusion, we have developed an efficient FEM- 
enabled numerical method for predicting the support- 
induced dissipation in micro- and nanoscale mechanical 
resonators. In combination with existing models fo r othe r 
relevant damping channels (e.g. fluidic and TEEPSlSlj^ 
our "phonon-tunneling" solver makes further strides to- 
wards accurate prediction of the mechanical quality Q. 
Furthermore, we demonstrate a stringent experimental 
test of the corresponding theory using resonators engi- 
neered to isolate support-induced losses. Our results 
demonstrate unambiguously that phonon-tunneling plays 
a significant role in the mechanical dissipation of these 
devices and illustrates the strong geometric character of 
this fundamental damping mechanism. Finally, we note 
that the weak coupling approximation underlying our 
treatment is more general than the condition of small 
contact area. Thus, our numerical solver can in principle 
be extended to other relevant scenarios such as systems in 
which the stresses at the contact are strongly sup pressed, 
as for example phononic-band-gap structures^. 
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I. METHODS 

Epitaxial Materials Structure and Resonator 
Fabrication Procedure - The materials structure for 
our high refiectivity resonators consists of 40.5-periods 
of alternating quarter- wave GaAs (high index) and AlAs 
(low index) layers grown lattice-matched to an off-cut 
monocrystalline germanium substrate. The ideal total 
thickness of the heterostructure is 6857.6 nm, with indi- 
vidual layer thicknesses of 77.6 nm and 91.9 nm for the 
GaAs and AlAs respectively; yielding a nominal peak 
reflectivity at 1064 nm, as with our previous optome- 
chanics experiments^. With this design, the germanium 
substrate enables the use of a high-selectivity gas-phase 
etching procedure, based on the noble-gas halide XeF2, 
in order to rapidly and selectively undercut the underly- 
ing germanium substrate. Thus we realize a free-standing 
epitaxial Bragg mirror via a simple and fast-turnaround 
fabrication procedure. The details of both the epitax- 
ial materials design and microfabrication procedure are 
covered in Ref. 48, 

Measurement technique - To characterize the fre- 



quency response of our microresonators we utilize a 
custom-built optical fiber interferometer featuring a con- 
tinuous fiow ^He cryostat as the sample chambeiHSl. 
High sensitivity displacement resolution is achieved in 
this system via optical homodyne interferometry. Cryo- 
genic testing of these devices is necessitated due to 
the limitations imposed by TED at room temperature. 
Estimations of the magnitude of TED is possible us- 
ing the analytical and finite element models developed 
previousljPSHEl^ which predict a Q-value of ~ 4000 for 
the current DBR composition and thickness at 1.8 MHz 
and 300 K — consistent with performed measurements. 
In order to avoid TED our cryostat enables interrogation 
down to 20 K (resulting in an estimated TED limited Q 
of 9.9 X 10^); the minimum temperature is currently lim- 
ited by the large view-port above the sample stage. Addi- 
tionally, this system is capable of vacuum levels down to 
2.5 X 10~ millibar at cryogenic temperatures, removing 
any additional damping induced by fluidic, or squeeze 
film effect82^"2Sl. The eigenmodes of the resonator are 
excited by driving a high-frequency (10 MHz) piezo disc 
soldered to a copper stage in thermal contact with the 
cold finger. For spectral characterization, the piezo disc 
is driven with white noise and the resonator frequency 
response is recorded on a low-noise spectrum analyzer. 
For the free-ringdown measurements, the decay of a res- 
onantly excited device is recorded in a single shot on a 
high-speed oscilloscope. 



II. SUPPLEMENTARY INFORMATION 

Appendix A: Numerical calculation of Q-values 

A rigorous derivation of Eq. (1) is given in Ref. [21] Al- 
ternatively, if one uses a decomposition of the displace- 
ment field in terms of the unperturbed support modes 
and the discrete modes of the resonator volume, Eq. (fl]) 
follows simply from applying Fermi's Golden rule to 
phonon decay with the interaction Hamiltonian between 
the resonator volume (labeled <) and the surrounding 
supports (labeled >) given by /^ dS ■ cr^ ■ u^ for case (i) 
and — JgdS ■ cr^ ■ Uy for case (ii). Within this frame- 
work, it is straightforward to realize that the validity of 
Eq. (Ill) is more general than the condition kud <C 1 and 
will also apply to systems in which the stresses at the 
contact are exponentially suppressed, like for example 
phononic-band-gap structures. In general the decompo- 
sition between "resonator volume" and "supports" con- 
sistent with the weak coupling condition need not be 
unique. For our case the use of this master formula is 
completely equivalent to prior intuitive approaches based 
on forcing the substrate with the stress source generated 
by the resonator, as can be shown rigorously by using 
— for the elastic Green's function of the substrate — a 
spectral decomposition in terms of its free modes. Fi- 
nally, one should note that Eq. (fTl) assumes that the reso- 
nance is spectrally resolved and that the support-induced 



frequency shifts are small compared with the free spec- 
tral range^I]. Thus, in the presence of mode-couplin^22|42| 
our treatment remains valid provided that the mode- 
coupling is not dominated by support-induced interac- 
tions, which includes the case where it is accounted for 
by FEM assuming perfect clamping and excludes cases 
where symmetry-breaking induced by the support is rel- 
evant. In the particular case studied here, the suspended 
structure has reflection symmetries with respect to the 
X — y, X — z, and y — z planes; while the entire structure 
including the substrate is only invariant under the latter 
two, allowing one to classify the modes accordingly. 

To derive Eq. ^ from Eq. ([I]) we adopt for the free 
elastic half-space, modelling the decoupled support, a 
decomposition into eigenmodes "),§„_, {r) (with n = 
0, ±1, ±2, . . . and g > 0) that have axial symmetry with 
respect to z (cf. Fig. 1). These are related to the plane 
wave eigenmodes Uqli{f) by 



,(0) 



',71,7 



(f) 



^/2^i J—^ 



d^(- 



-i)" e™^fi(°' 



'lilfi^v),! 
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where j = I, t, s labels the different types of rele- 
vant modes [i.e. longitudinal (/), transverse SV (i), and 
SAW (s) given that SH waves do not contribute] and 
we adopt spherical coordinates for the incident wavevec- 
tor q{q,9,ip) = q{sm 9 cos ip , sin 9 sin ^ , cos 9) [9 = 7r/2 
for 7 = s and 9 < 7r/2 otherwise]. We note that for 
the suspended plate geometry considered, the appropri- 
ate resonator mode satisfies m'^ = at the contact S so 
that we need to evaluate the second term in Eq. ([l]). 

The thin plate condition t -^ R directly allows 
us, given the flexural nature of the modes of interest, 
to neglect stresses at S that are parallel to the sub- 
strate with the possible exception of bending-momcnt 



contributiongSS — this also applies if there are small 
transverse dimensions comparable to t. However the 
bending-moment contributions also become negligible in 
the limit t/R — > as can be shown by using: (i) that 

given uj]i <C c^/R V7, we can Taylor expand u„ g ^ (f) 
at the origin in the integral over S, (ii) that we can as- 
sume relevant stresses to be concentrated around the 
ends of the auxiliary beams so that the bending mo- 
ments at S are mostly oriented along y, (iii) the reflection 
symmetries, and (iv) that, barring interference effects, 
these ben ding -moment contributions are at most of rel- 
ative orde^^knt — here kn = {I2pii/Euy/^ y^uju/t is 
the resonant wavevector for the propagation of flexural 
waves. Thus, we find that for all mode types other than 
— h [antisymmetric (symmetric) with respect to Rx (Ry)] 
the correction associated to neglecting the bending mo- 
ments scales as AQ/Q ^ (fc^^t)^, while for — h modes 
it scales as AQ/Q ~ /c/jt (note that L ^ R). In turn, 
we find that the relative error in using Eq. (fTl), arising 
from the weak coupling approximation, scales in this case 
as AQ/Q ^ \Aui^\/ujR -- \Ai{ijJr)\/ur ^ {kntf., where 
the phonon-tunneling induced frequency shift is approxi- 
mated by Ai{LLiii) K, —{2/tt) L AijjI{ijj)/ijj with I{lu) the 
environmental spectrum^^. 

Hence we can assume dS ■ cr'j^ \\ z and neglect the 
variation of fr I (f) across the thickness t (i.e. the 
z-dependence at S), so that the support modes only 

enter into Eq. M through u]^ g ^ {f)\z=o- To deter- 
mine the latter we adopt cylindrical coordinates f = 
{r cos (j),r sin (t>,z), substitute into Eq. (Al) the plane 



wave expressions for u-J g ^ ■zi''")^ exploit that reflec- 
tion at the free surface preserves the tangential com- 
ponent of the wavevector, and use the Besscl integral 
Jnix) = ^/_^^e-*("'^-^'''"'^)d0. Thus we obtain 
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that substituted into z ■ Uq[r) in Eq. (1) leads to Eq. (2) after using J — > X^n 7 /o°° ^9?'' 

C d9{T: /2)'^i^^{sin9Y'<^'^9 [d^ is the dimensionality, i.e. d^ = 2 for 7 = 5 and ^..^ = 3 for 7 7^ s), performing 
the substitution uj = c^q (for each 7), and integrating over w. Subsequently, substitution of the explicit expressions 

for the plane wave eigenmodes Ug' (f) (see for example Appendix A in Ref. I2ip and v — cos 9 into the deflnition of 
Uji,y{q,Vs) allows us to obtain: 
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where we use the ratio a = (ctjci)^ = {\ — 2vs)/2{l — Vs) 
for the supports' material {vg is the corresponding Pois- 
son ratio). In turn ^(a) is the ratio of the propagation 
velocity of surface waves to C(, which is a function of a 
that is always less than unity^, and 



C{a) 



2-aC\a) 



[e{a)/4 + e(a)~l] 



^1/2 



(A4) 



Vl - aeia) j 
The sum in Eq. (pi) can be reduced to a sum over n > 



by noting that J_„(a:) = {—l)Jn{x) and that as the res- 
onator mode is real, the linear stress Fourier components 
satisfy /z,-ti — It n- Furthermore, the length of the cen- 
tral resonator L is comparable to the radius R and we 
focus on low-lying resonances of the suspended structure 
so that the aforementioned condition ur <S^ Cj/RVj is 
always satisfied. This implies U2'y'^rn,j/J2^''^n,j\ ^ 1 
for m > n and ^ 



'.y '^n,7 



7^ 0, which can be understood by 
considering the behavior of the Bessel functions for small 
arguments. Thus we find that in Eq. (pi) the sum over 
the index n is dominated by the first non- vanishing term 
as determined by the reflection symmetries Rx, Ry The 
latter also imply (n = 0, 1, 2, . . .): 



f+2n+i =R j d(t) Idzf- (T+f (f) • z sin(2n -|- 1)0 , 
fz.2n+i =R j d(t) dzf- (T]^*{f) ■ z coa{2n + l)(j> , 



fz,2n=R d(j) dz f • (T^ (r) • z sin 2n</) , 



Jz,2n+1 ~Jz,2n ~ Jz.2n ~ Jz,2n^ 



0; 



(A5) 



where the resonator mode of type a, P satisfies Rxcr'^ — ctcr'^ and Rycr'^ — (3(t'^ . To efficiently extract the 
above from the FEM simulation we convert them into volume integrals using an adequate Gaussian weight so that 
for example, for a fully symmetric mode we have 



fz.L = lim. 



a,-^o V7ra 



* Jv 



dr^e V "* '' f ■ cr'j^^r) ■ z cos 2n(j) , 



(A6) 



where we again use cylindrical coordinates and V denotes 
the resonator volume. In addition we exploit that the re- 
flection symmetries naturally allow to perform the FEM 
simulation on a single quadrant. Thus, numerical eval- 
uation can be conveniently performed using a fixed a* 
and a mesh size M such that (V/AM)''-/^ < a* < t. We 
have checked the convergence and estimate the numerical 
error to be of order 5%. 

Numerical simulations of the resonator mode are 



performed with the aid of COMSOL multiphysics. 
Accurate three-dimensional CAD models representing 
the resonator geometry are generated using Solidworks 
(matched with high quality SEM images as described in 
Appendix^ and the bidirectional interface between the 
two programs is exploited in order to perform a paramet- 
ric sweep of the auxiliary beam contact position for de- 
termining the pertinent information about the relevant 
mode; namely, its eigenfrequency, linear stress Fourier 
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components fz,n, and normalization constant. In this in- 
stance a single CAD file is used with a global variable 
incorporated in order to control the lateral position of 
the auxiliary beams with respect to the center of the 
central resonator. We use for the mechanical properties 
of our single-crystal resonators an anisotropic materials 
model incorporating the elastic stiffness matrix for the 
epitaxial structure as obtained from a weighted average 
between the relative content of GaAs and AlAs (46.37% 
GaAs / 53.63% AlAs). The corresponding parameters 
are: Cn = 119.6 GPa, C12 = 55.5 GPa, C44 = 59.1 GPa, 
and PR = 4483 kg/m'^. The resonator axes are aligned 
along (100) (zinc-blende structure). Note that we ignore 
the 6-degree misorientation of the germanium substrate 
as we have checked that it has a negligible impact (er- 
ror of 0.3%) on the simulated frequency response of the 
resonators. 

We note that our generic symmetric inscribed structure 
includes the particular case of bridge geometries with no 
undercut for which a simple variant of the method used in 
Ref. |2T] allows us to obtain an analytical approximation 
for the Q-value of the fundamental mode Qc~c valid for 
37rt/2L <^ 1. Our scenario differs from the one considered 
in Ref.[21jin two ways: (i) there is now a single half-space 
support instead of two and (ii) its free surface is oriented 
parallel to the beam's axis instead of perpendicular to it. 
Thus, for the fundamental flexural mode as the resonant 
wavelength in the support is much larger than the bridge 
length L, to lowest order in (t/L)^ the stresses at both 
clamping points add coherently so that the overall effect 
of (i) is to double the dissipation. In turn (ii) implies that 
the roles of the dimensionless displacements for compres- 
sion and bending are interchanged so that the dissipation 
of the vertical bending modes are further corrected by a 
factor oiuc/uv Thus if we consider that the support and 
resonator are made of the same material characterized by 
a Poisson ratio i^ = 1/3 we obtain 



^c—c 



0.92 L^ 

TT^Wt* 



(A7) 



Hence as a non-trivial check we have applied our nu- 
merical method to the fundamental mode of clamped- 
clamped square beam monolithic geometries with no un- 
dercut, u = 1/3, and aspect ratios L/t ranging from 15 
to 40 and compared the results with those correspond- 
ing to Eq. (A7). We find a discrepancy e that decreases 



monotonously from 20% to 4% which is consistent with 
the rough heuristic estimate e ^ 37rt/2i. 



Appendix B: Analysis of Completed Devices 

The resonator layout we have designed features 16 de- 
vices on chip, each with identical central resonator di- 
mensions (nominally 130 x 40 /im^). The 16 devices are 
divided into two sub-units featuring different outer radii 
(116 um and 131 um respectively), which are included in 
order to probe the effects of the auxiliary beam length on 



the dissipation. Finally, each of the two subsets contains 
8 variations of the auxiliary beam contact position, vary- 
ing from the center to the extreme outer edge of the cen- 
tral beam, with a single design chosen to match the the- 
oretically calculated node position (auxiliary beam posi- 
tions of 13, 21, 29, 37.4, 44, 50, 56, and finally, 62.5 ^m). 
To ensure a thorough investigation of each geometry, two 
separate but nominally identical chips are measured. 

Central to this study is an accurate determination 
of the geometric properties of the optomechanical res- 
onators. Thus, we employ a variety of analytical tech- 
niques for the characterization of these devices as detailed 
below. We find that the actual thickness of the DBR is 
6.67 /zm, the central resonator dimensions are enlarged 
by 1 /xm at each free edge as compared with the nominal 
design values, and finally, Lund is destructively measured 
post characterization and found to be on average 27 fim. 
In turn the microfabrication procedure detailed Ref. |48| 
entails h ~ iund ■ 

Thickness: In order to accurately determine the physi- 
cal thickness of the resonators, we rely on measurements 
of the DBR reflectance spectrum. This procedure be- 
gins by recording the reflectance of the mirror stack (on 
wafer) as a function of wavelength via spectrophotom- 
etry. A transmission matrix model is then used to fit 
the measured high-reflectivity stop-band; the individual 
layer thicknesses are adjusted assuming constant (fixed 
percentage) growth errors for the constituent films. Note 
that the wavelength of peak reflectivity of the mirror is 
highly sensitive to variations in layer thickness. In fact 
for this structure, a 1 nm variation in the individual layer 
thickness shifts the wavelength of peak reflectivity by ap- 
proximately 10 nm. Relying on accurate knowledge of the 
room temperature refractive index of the binary films, we 
realize a minimum wavelength resolution of ±lnm; thus, 
the thickness accuracy is better than 20 nm for the DBR. 
From this analysis we have determined that the actual 
thickness of the DBR is slightly shorter than desired at 
6.67/im with a peak on- wafer reflectance near 1060 nm at 
room temperature (ideal target thickness of 6.86/xm, cor- 
responding to a peak wavelength of 1078 nm at 300 K). 
The thickness is further verified by scanning probe mea- 
surement of the DBR following the anisotropic etch of 
the epitaxial layers. The profilometer provides an upper 
limit to the DBR thickness, as additional etching arising 
from surface sputtering of the Ge substrate is unavoid- 
able. These measurements yield a conservative thickness 
estimate between 6.7/im and 6.8/im, verifying the more 
accurate spectrophotometer derived value. 

Resonator Dimensions: The lateral dimensions of the 
resonators are determined by obtaining high resolution 
micrographs of each individual structure in a field emis- 
sion scanning electron microscope (Zeiss Gemini). Im- 
age analysis shows that the lateral dimensions of the 
resonators have expanded by -1-1 ^m on each edge, with 
the following results: reducing the nominal external sup- 
port diameter by 2 /zm, increasing the auxiliary beam 
width from 5/im to 7/im, and increasing the overall 
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lateral dimensions of the resonator by 2 fim to 132 /im 
and 42 /xm. Additionally, a combination of process non- 
idealities (non-optimized exposure or development times) 
during lithography result in the formation of a 3-/im- 
radius fillet of at each corner of the device. These results 
are fed back into the CAD model of the resonator in order 
to generate the true resonator geometry for simulation. 
An overlay of the simulated resonator geometry and mi- 
crographs obtained via scanning electron microscopy can 
be seen in Fig. 3(a). Note that the resonators used in this 
study were not subject to potentially damaging energetic 
processes beyond the required plasma etching, including 
both SEM and FIB (as described below), until all dissi- 
pation measurements had been completed. 

Undercut: In order to perform measurements of the 
support undercut distance, a dual beam SEM/FIB (Zeiss 
Gemini) is utilized to mill a window through the DBR 
and expose the underlying germanium. Because the 
GaAs/AlAs heterostructure is opaque to visible light, 
it is not possible to simply view the undercut distance 
with an optical microscope. This method allows for an 
accurate determination of the lateral etch distance be- 
low the supports. Image analysis yields an average dis- 
tance of 27 fiUi for the structures. Note that multiple 
chips of identical geometry were released simultaneously 
in a single process run, in order to ensure repeatability in 
the resonator dimensions. Selected measurements across 
the chip verify that the undercut length is constant for 
the resonators studied here (measured values fall between 
26.5 and 28.2 fim). 



Appendix C: Q- value and frequency measurements 

We utilize two options for driving and characterizing 
the resonance of interest: (i) by applying broadband 
white noise to the piezo disc for extraction of the me- 
chanical frequency spectrum (simultaneously driving all 
modes within the system bandwidth), and (ii) by excit- 
ing a desired mode resonantly with a sinusoidal voltage 
input, abruptly shutting off the drive, and then recording 
the free-ringdown of the structure. In the first method, 
Q is extracted by measuring the width of the resonance 
of interest, while in the latter, the single-shot amplitude 
decay time of the ringing structure provides the damping 
rate of the resonator. 

For data analysis, we employ a combination of spectral 
fitting with a Lorentzian function, with linear-regression- 
based fitting of a decaying exponential in the case of the 
ringdown data. The envelope of the raw ringdown sig- 
nal is created by first squaring the dataset (in order to 
utilize both the positive and negative components of the 
decaying sinusoid) and then averaging over a 10-20 period 
window (it is important to note that a typical ringdown 
dataset contains more than 1.5 x 10'^ periods of oscilla- 
tion). Finally, to linearize the data, we simply take the 
natural logarithm of the mean-squared amplitude. In 
contrast to the single-shot ringdown datasets, the spec- 
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FIG. 4: Measured bistable response of the anti-symmetric 
mode. 



tral measurements require multiple averages for a clean 
signal (typically ^ 50 — 100). The Lorentzian fit param- 
eters include the amplitude, center frequency, and full 
width at half maximum (FWHM), with the latter two 
used for calculating Q from their ratio. The fast Fourier 
transform (FFT) of the ringdown signal can also be used 
to determine the eigenfrequency, while the 1/e decay time 
T allows for the extraction of the resonator quality factor 
via the relation r — Q/nf. 

Mode identification is realized by comparing the res- 
onator frequency response as a function of geometry with 
the simulated eigenfrequency values. The modes are fur- 
ther distinguished by the relative geometric-induced non- 
linearity at resonance. The desired free-free mode re- 
mains linear to the limits of our piezoelectric-based in- 
ertial drive. On the other hand, the neighboring anti- 
symmetric mode exhibits a significant hardening spring 
Duffing response and can readily be driven into a bistable 
regime as shown in Fig. |4j Care is taken to drive this 
mode below the threshold for bistability to avoid compli- 
cations in dissipation extraction. This marked difference 
in the responses of these two types of modes is consistent 
with their free-free versus clamped-clamped nature^. 

To each dissipation mechanism there is an associated 
dispersive effect induced by the interactions with the cor- 
responding environment that shifts the resonant frequen- 
cies. For a given resonance {ujji) this shift can be posi- 
tive or negative depending on whether the environmental 
spectrum is dominated, respectively, by modes with fre- 
quencies lower or higher than ujj^. In turn, the two types 
of modes have markedly different surface-to- volume ra- 
tios, larger for the antisymmetric resonance and smaller 
for the symmetric one; exhibit a positive shift which is 
substantially larger for the symmetric mode, and a "back- 
ground" dissipation that is larger for the antisymmetric 
one (cf. Fig. 3). These facts can be reconciled by assum- 
ing that two materials-related dissipation mechanisms 
contribute to the "background": a bulk one leading to 
an overall positive shift which is the same for both types 
of modes, and a surface one leading to a smaller negative 
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shift that naturally scales as the surface-to- volume ratio. 
In turn, theoretical estimates for the phonon-tunneling 
induced shift yield a negligible negative shift that should 
also follow the mode profile leading to a significant mod- 



ulation that is not observed. Likewise mode coupling be- 
tween the different resonators would also be incompatible 
with a constant shift. 
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